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ABSTRACT 

The principle of multiple material layering is well known as an effective 
method of reducing heat transfer; however, thermal gradients can impose sig- 
nificant mechanical loads and lead to delamination and subsequent component 
failure. An analysis is developed and experimental data are discussed for 
the thermomechanical effects of multilayered materials on a heat sink sub- 
strate of cylindrical geometry subject to thermal cycling. The geometry is 
heated in cross-flow by a high-velocity flame and cooled in crossflow by 
ambient-temperature air from a critical flow orifice. Each layer of material 
possesses a threshold beyond which small changes in temperature or mechanical 
loading greatly influence the life in thermal cycling of-the layered mat- 
erials. Comparisons are made between the thermomechanical loads predicted by 
various numerical codes for the linear case and by a simplified analytic 
model. 

NOMENCLATURE 

A area 

AR aspect ratio 

a,b exponent, eq. (8) 

Co constant, eq. (8) 

D diameter 

E modulus of elasticity 

F temperature ratio function, eq. (2) 

FE finite element 

h heat transfer coefficient 

k thermal conductivity 

L axial length 

M Mach number 

N number of cycles 

Pr Prandtl number 

q heat flux 

Re Reynolds number 

r radial coordinate, rads 

So standard deviation 

T temperature 

t thickness 



u displacement 

y exponent, eq. (2) 

z axial coordinate 

a thermal expansion coefficient 

a differential 

e strain 

e strain rate 

e plastic strain 

0 " circumferential coordinate 

0 stress 

t time 

Subscripts : 

cr creep 

h heating 

f fi lm 

p pulloff 

q heat flux 

s surface 

rr radial 

zz axial 

zr axial-radial (shear) 

ee circumferential 

® free stream, bulk 

0 reference value 

1,2 reference materials 

Superscripts: 

+ transient, eq. (19) 

* related to pull-off, eq. (19) 

INTRODUCTION 

Toward Higher Engine Efficiencies 

Improved engine efficiency without a sacrifice in operating life is a 
major goal of the propulsion community. In 1981, U.S. airlines consumed a 
billion gallons of fuel; a one-point increase in efficiency could have re- 
sulted in a savings of some $60 M (ref. 1). One method of achieving higher 
efficiencies is to increase the cycle operating temperature. To prolong 
operating life at these temperatures, components can be covered with multiple 
layers of plasma-sprayed materials. 

Plasma Sprayed Surfaces 

The potential use of plasma-sprayed layered materials is promising, but 
not without difficulties and limitations. Previous experimental work 
(refs. 2-8) has shown that yttria-stabilized zirconia (Zr0 2 -Y 2 C>3) ceramic 
coatings plasma sprayed onto a metal substrate fail after varying lengths of 
time when cyclically heated by exposure to high-temperature, high-heat-flux 
flames. The ability to predict the life of these layered materials is made 
difficult because several modes of failure or degradation exist, including 
spalling, erosion, corrosive attack, and corrosive oxidation of the bondcoat. 
The failure occurs principally in the ceramic at the bondcoat interface. 
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Although the mechanisms that control these failure modes are not well under- 
stood, when spalling or debonding occurs, the adverse effects include direct 
exposure of structural substrates to the high-temperature environment. 

Objectives 

In many cases the heat transfer can be evaluated, but the coupled ther- 
momechanical response of specimens in both the spacial coordinates and time 
is not known. One needs to investigate the effects of heating and cooling 
transients, time at temperature, and heat flux, coupled with elastic and 
inelastic material behavior. 

In this paper we model the cyclic heating and cooling of a cylindrical, 
multilayered specimen alternately placed in high-velocity hot and cold gas 
streams, in crossflow. We want to address variations in thermal profiles 
and stresses due to modeling and computational codes by having researchers 
independently, but cooperatively, look at the same problem. Thus we look at 
the results from the thermal codes ADINAT, MARC, and SINDA and the stress 
codes ADINA, MARC, and FEATS. 

MODELING CONSIDERATIONS 

Thermocycling 

It has been determined that thermal stresses in layered protective 
ceramic coatings are sensitive to the heating rates and that short heating 
cycles greatly decrease life in terms of total time at temperature (ref. 2). 
Furthermore a 0.076-cm-thick surface material failed immediately and a 
0.038-cm-thick surface material failed some number of cycles later. Thus 
for our problem the heating and cooling boundary conditions will be related 
to short cycles, 3 to 5 minutes of heating, and 2 to 3 minutes of cooling. 

The thickness of the multilayered materials can be varied. For comparison, 
we start with a 0.038-cm thickness and thermocycle to an equilibrium temper- 
ature of 1311 K (1900° F). 

Limiting Life Temperature 

In reference 6 the effects of heating rate and operating temperature 
were explored by varying the initial or startup fuel-air ratio (F/A). At a 
reduced (F/A), flat-plate specimens ran over lxlCP cycles without failure, 
but at an elevated (F/A) such specimens ran less than 400 cycles to failure. 

Most importantly, the surface temperatures were 900 (1650) and 982° C (1800° F), 
respectively. The analysis showed a lower stress (20 percent) at 900° C, 
which did not explain the large change in cyclic life. This implies the 
existence of a threshold temperature for the model (i.e., above 900° C). 

Cycle Life and Adhesive/Cohesive Strength 

In reference 7 a statistical study of 22 cylindrical rod specimens (ex- 
posed to high-heat-flux conditions) and 24 cylinder disk pulloff specimens 
(no exposure) was made in an effort to define a failure model and relative 
strengths. Four failure mechanisms (strain relief, heating rate, stress 
reversal, and ceramic voids) were cited. The data sets were better described 
by log-normal distributions. Such a distribution where the failure rate 
decreases as run cycles increase could be related to the threshold tempera- 
ture, as stated in the preceding section. The data of reference 7 will be 
used in our calculations. 
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Effects of Combustion Gases 

In reference 8 tests with salts indicated that some would cause immediate 
coating failures (e.g., vanadium compounds), while others would take somewhat 
longer (e.g., magnesium compounds). The failure mechanism could be the in- 
terstitial condensation of salts in the ceramic along with oxidation of the 
bondcoat. In this analysis, we consider only clean gases. 

Substrate Precooling 

In reference 9 the substrate temperature was systematically varied, and 
in reference 10 an optimum region 21 to 316° C was reported. This implies 
some residual stresses, but their magnitudes are unknown, so we will consider 
an initially stress-free coating. 

PI asma Arc Power 

In reference 11 the use of plasma-spray arc power increased specimen 
life somewhat with increasing arc current (525 to 950 A). At least two ef- 
fects are noted; the porosity of the coating decreases and the development 
of a microcrack network increases with increasing arc power (See also the 
motion picture supplement). In addition, a higher quench heat transfer 
coefficient due to increased impact velocity, substrate heatup, larger 
droplet mobility, and wetting characteristics upon impact are suspected to 
have an effect. We will not consider microstructure in our model directly; 
however, we will use the sample data sets of reference 7, which were sprayed 
at 650 A. 

Thermal Shock Effects 

Crack propagation and arrest (thermal shock) are inherent to all these 
studies. Analysis and testing of Si O 2 and Zr02 fiber materials demon- 
strated high-temperature capabilities, but strengths are low and bonding 
techniques undeveloped (ref. 12). 

Since no microstructure will be considered, the layers modeled depend 
only on bulk material properties. However, in this analysis we will consider 
only one type of coating, Zr 02 -V 203 (yttria-stabi lized zirconia, or YSZ). 

Inelastic Behavior (Creep) - A Significant Factor 

Some effects of high-temperature inelastic behavior are presented in 
reference 13, where a form for creep behavior is established. Furnace tests 
demonstrate creep effects for plasma-sprayed YSZ disk specimens and YSZ-Bond- 
Substrate bar specimens. These tests show that significant creep begins near 
1000° C (see the section Limiting Life Temperature) and that the bondcoat 
(NiCrAlY) still retains sufficient strength to creep the YSZ at temperatures 
to 1177° C. 

In all cases the experiments and analysis show a sensitivity to tempera- 
ture and thermophysical properties. These results may be qualitatively sum- 
marized as: ' 

1. That if the sample can withstand the thermal shock (remain within the 
elastic range), then creep and oxidation/corrosion must be considered. 

2. At increased thermomechanical loading, failure can occur during cool- 
down or subsequent thermal shock phases where the additional strain due to 
creep causes a stress load beyond the capacity of the coating and it fails. 

For the initial development of the numerical models, elastic behavior is 
considered; in the simplified model, creep is considered. 
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Toward Ceramic Components 

A monolithic ceramic turbine blade with cooling slots was proposed, and 
operated, during the World War II; see reference 14 for photographs. More 
recently, a similar concept was proposed in reference 15 using YSZ. Compu- 
tational heat transfer and industrial applications of ceramic-sprayed com- 
ponents (turbine blades, seals, pistons, etc.) are explored in references 16 
to 19. The details of predicting turbine blade temperatures and the narrow 
operating range (within ± 30° C) along with the thermal profiles for the 
sprayed blades are provided in references 20 and 21. In reference 22 various 
heat transfer correlations for internal cooling of turbine blades are inves- 
tigated, and recommendations and limitations for their use are presented. 

Although the model can be applied to these more complex configurations, 
we will model a cylindrical geometry. 

HEAT TRANSFER ANALYSIS 


Heat Transfer Coefficients 

Accurate information for the convective heat transfer coefficient to and 
from cylinders exposed in a stream of air is essential for evaluating the 
cylinder's surface temperature and heat flux to and from this surface. In 
addition to the convective path there is a radiative path, but the latter 
becomes important only at relatively low Reynolds numbers (lxlCr to lxlCP) 
or when rather large temperature differences are acting between the cylin- 
der's surface and the environment. 

From a literature study (refs. 23-27) we have decided to adopt separate 
heat transfer coefficients for the heating and cooling of the cylinder. The 
correlations are valid for cylinders with smooth surfaces, and will be used 
even though the plasma-sprayed finish of the ceramic cannot be considered 
smooth. 

Heating . A relationship put forward by Churchill and Brier (ref. 23) 
which correlates the experimental data (refs. 27-29) with the physical prop- 
erties evaluated at T* was adopted. 

Although the test cylinder temperature varies in both the axial and cir- 
cumferential coordinates, the variations are small relative to the radial 
gradient. Thus, integrating over e, the average value of the Reynolds 
number exponent becomes n = 1/2; the model is now assumed to be axi sym- 
metric . 


Nu = Ip- = 0.6 Re 1/2 Pr 1/3 

CD K CD OO 



( 1 ) 


According to Morgan (ref. 27) when the intensity of the turbulence is 
taken into consideration for Re > 1x10^, the heat transfer coefficient 
can increase by 18 to 37 percent depending on the intensity. It is apparent 
that neither reference 23 or King (ref. 29) nor Hiplert (ref. 30) have in- 
cluded the effects of the turbulence intensity in their correlations, and 
this also represents a limitation of our analysis. 


Cooling . For the cooling part of the cycle, most of the experimental 
data available have been fitted to a family of equations over a range of 
Reynolds numbers. Fand and Keswani (ref. 24) developed correlation for a 
Reynolds number range between lxlO - ^ and 2x10^. The expression they 
proposed is 
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( 2 ) 


Nu f = ( 0.184 


+ 0.324 Re* /2 + 0.291 Re°* 5+y )l 


/ T f\ 
= \'T 

\ no ' 


0.17 


F 

where 


y = -0.253 + 0.0407 Re 8 ' 168 


Tf and T® being the absolute film and ambient temperatures, respectively. 
Since the properties are evaluated at the film temperature, the heat trans- 
fer coefficient is no longer represented by a simple algebraic function in 
terms Ts- 

The third term inside the brackets of equation (2) represents the heat 
transfer on the front and rear portions of the cylinder; and equation (2) 
calculates an averaged heat transfer coefficient for the whole cylinder. 

NUMERICAL THERMAL MODELING CODES 

In general, the experimental applications will involve nonsymmetric geom- 
etries with temperature gradients in all three coordinates (r,e,z) due to 
external and internal flows. Herein we will treat the body as a uniformly 
heated cylindrical heat sink. 

ADINAT (ref. 30) . 

A 0.038-cm-thick ceramic coating, yttriastabilized zirconia (YSZ), is 
plasma sprayed over a 0.013-cm-thick nickel-chromium-aluminum-yttria 
(NiCrAlY) bondcoat. The stainless substrate* is 1.30 cm in diameter.** 

Figure 1(a) shows the finite element model of the system of 18 elements 
and 93 nodes. Although the 4-node element could be used, the 8-node element 
was chosen in order to obtain better results for the stress calculation. A 
refined grid with 333 nodes and 66 elements produced nearly the same temper- 
ature distribution. 

Heating and Cooling Cycles . Figure 2 presents temperature distributions 
during the heating and cooling periods. 


*The substrate material used in the tests (e.g., refs. 2, 22) Rene 41 or 
inconel are classified as super alloys. Upon direct exposure to the 0.3 
Mach burner flame, these substrates corrode slowly; 304-stainless steel cor- 
rodes rapidly, has a higher thermal expansion coefficient and would be a 
poor choice of substrate material; the coating is, however, very protective 
and 304 is used in this study for simplicity and as a parametric term. 

**1.27 cm diameter was used in the MARC code and both 1.27 and 1.3 cm were 
used in SINDA-FEATS with little change in the results. 
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The hot gas temperature is 1311 K during the first 5 minutes, and the 
cold gas temperature is 298 K in the next 3 minutes. The label "Zr02- 
interface" represents the ceramic-bondcoat interface. Since the bondcoat- 
substrate interface temperature difference is only about 5 K, it is not 
shown. In the calculations, all thermal properties except density are con- 
sidered to be temperature dependent. Those properties, shown in table I, 
are the input to ADINAT. The program picks up the appropriate value by 
linear interpolation. 

For heating, the heat transfer coefficient proposed by Church and Brier 
(ref. 23, eq. (1)) was used. Since all properties are evaluated at the bulk 
stream temperature of the hot gas, equation (1) can be easily converted to a 
simple formula with the aid of gas tables and known burner Mach number (0.3); 
for our problem this becomes 


916.6 

" t 0.12 


W 

9 

m-K 



This formula is also fed into ADINAT through input. 

As for the cooling process, the correlation by Fand and Keswani (ref. 24, 
eq. (2)) has been adopted. All fluid properties are now evaluated at the 
film temperature, and the heat transfer coefficient is a more complex func- 
tion of surface temperature. A table of h vs T s was constructed and 
was also fed into ADINAT through input. 


Coupling of the Heat Transfer Coefficient . To demonstrate the significance 
of a variable heat transfer coefficient, we repeated the same problem of 
heating for the first 14 seconds by using a uniform heat transfer coeffi- 
cient, h = 511 W/m^-K (90 B/hr-ft^-F). The dashed lines in figure 3 
represent the numerical solutions for this case. As can be seen, the dif- 
ference between the solid and dashed lines continues to grow as time in- 
creases, until equilibrium is reached. Since the temperature of the solid 
varies over such a wide range, and in turn gives rise to thermal stress 
fluctuations, the temperature-dependent heat transfer coefficient and other 
thermal properties should not be neglected. 


MARC (ref. 31). 

For the axisymmetric heat transfer analysis, internal to MARC, the geom- 
etric configuration discussed in connection with ADINAT was adopted. This 
method, like ADINAT, is a finite element method. 

Figure 1(b) shows the finite element model consisting of 72 elements and 
293 nodes. The 8-node isoparametric, axisymmetric element was used to obtain 
both thermal and stress calculations. The system simulated a set of 4 
stacked disks (elements) when combined with a symmetry condition along the 
center plane; this approximates a cylinder of length 2L. 

Although the thermal analysis is essentially one-dimensional, the use of 
several elements in the axial direction allows coupling of the thermal and 
stress analyses without the use of a mapping function from one mesh to 
another. 

In the section Temperature profiles , we cite the use of symmetry and 
free end conditions on z to simulate a short cylinder and fixed end con- 
ditions at both z = 0 and z = 2L to simulate the infinite cylinder. For 
either case the temperature does not vary in the axial (z) direction. 

Thermal boundary condition subroutine . A subroutine was written to 
describe the convective (radiative) boundary conditions. These conditions 
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have been described above and are given by equations (1) and (2) for heating 
and cooling, respectively. Either the Dirichlet or Neumann conditions can 
be calculated in such a subroutine as a normal part of the MARC thermal 
analysis and can be easily adapted to handle more complex formulations such 
as may be required for axial and circumferential variations in surface tem- 
perature. 

Temperature profiles . Graphs of the surface, the interfaces between the 
ceramic, bondcoat, and substrate and the centerline for heating and cooling, 
are similar to those shown as figure 2 and will not be repeated here. The 
difference here is that we have elected to use a superalloy (ref. 17) to 
represent the substrate. Since the superalloy has a lower thermal conduc- 
tivity than type 304 stainless steel, the relative interface temperatures 
are about 50 K higher for the transient. 

SINDA (refs. 32,33). 

In contrast to the other codes, which are finite element, SINDA is a 
finite difference code. Again, the geometry adopted is that described in 
the section Heat Transfer Coefficients. The properties are now represented 
as lumped values. To make the temperature profiles available to the mechan- 
ical code FEATS, the average element temperature must be interpolated to 
provide values at each node. Across each interface the average thermal con- 
duction is represented by the inverse sum of the conductivities. 

We used both a 40-element, 82-node model and a 105-element, 212-node 
model and achieved about the same results with both. The models, shown in 
figure 1(c), had 10 elements in the ceramic, 10 in the bondcoat, and 10 in 
the first 0.038 cm (0.015 in.) of the substrate. The remainder of the ele- 
ments, either 10 or 75, uniformly partitioned the remainder of the substrate. 
It was felt that such a division would capture the extreme gradients in the 
ceramic, especially at the bondcoat-ceramic interface. 

Comparisons 

Although three different numerical techniques, average and variable 
thermophysical properties, and heat transfer coefficients were explored, the 
profiles were within 100 K. In this analysis, such a variation is small 
(10 percent) when compared with the absolute level and range of parameters 
considered. These results are also in reasonable agreement with the experi- 
mental temperatures measured at the interface (ref. 2). 

STRESS ANALYSIS 

The coupled heat transfer and stress analysis, termed thermomechanical 
analysis, enables one to combine the effects of mechanical loads imposed by 
surface tractions (e.g., skin friction and pressure drop) and body force 
effects (e.g., rotation, takeoff, and maneuvering) with those thermal loads 
imposed by changes in internal energy and volume. (It is hoped that the 
composition remains fixed.) Although these effects must be considered in 
application herein, we will consider a traction and body force-free axisym- 
metric geometry. 

Simplified Model 

Two-layer model with creep . Consider a two-material, cylindrical, heat 
sink model of multilayered ceramic with a stress-time history as illustrated 
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in figure 4. To determine the stress at the interface, including the effects 
of creep, the following equations must be satisfied: 

Compatibility - e? = 0 (4) 

Equilibrium =0 (5) 


Constitutive 


e = a AT + + 


f V 

J e cr 


For the simplified model (fig. 4(b)) the stress becomes 


( 6 ) 



Further assume that the system is isothermal (i.e., aT^ = AT 2 = aT) and 


„ a b -Q/RT 
e cr = Coo t e 


(8) 


Homogeneous body approximation . Let material 2 be the YSZ ceramic and 
material 1 the substrate, say Inconel. With the geometric and material con- 
stants given in the next section, A 2 /A 1 = 4t/d = 0.12, A 2 E 2 /A 1 E 1 = 0.07 


and assuming e cr 



>> e 


cr 


1 


equation (7) may be approximated as 


= | = -AaAT - 



(9) 


The ceramic creep rate may be expressed as (ref. 13). 


cr 


2.39xl0 _3 T _0,45 a 


-1230/(T-T ref ) 


( 10 ) 


where the units of x, a, and (T - T re f) are minutes, MPa, and K respectively. 
Differentiating, equation (9) becomes 


cr = -AqOT 

where 


-0.45 


-1230/ (T-T) 

A q = 2.39x10 J Ee e 


(ID 


( 12 ) 
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Solving, equation (11) becomes 


o 


e - E - 


e 0 


-1.82A 0 t 0,55 

e 


(13) 


From the initial conditions it is known that 

Gq = -a<xaT - Gp (14) 

and subtracting the initial strain equation (12) becomes 


£ 



-1.82A t 
e u 


0.55 


) 


(15) 


Effect of N cycles . For N heating cycles, each of time minutes, 
we may approximate the total time of creep as t = z-r-j ~ t^N. If the 
limiting strain is governed by a relation to the pull off strength of the 
ceramic (i.e., the adhesive-cohesive strength), the developed strain can be 
related as 


V 


(16) 


Estimates for interface stress . The following numerical values are used 
to estimate the stresses at the bondcoat-ceramic interface: 


Tc - To = 1000 K (1800° F) 

T s = 1311 K (1900 F) 

A« = 7.2xl0~ 6 (4xl0 -6 ) 

E = 1.38xl0 4 MPa (2x10° psi ) 

E 0 = 4.8xl0 4 MPa (7xl0 6 psi) 
r = 0.66 cm (0.26 in. ) 
t = OJ338 cm (0.015 in.) 
a p = op ± So (9 ±4.5 MPa) (1300 ± 650 psi) 

_tor 24-specimen set (ref. 7) 

°P = °D * ^0 

= 6.2 ± op/2 MPa (900 ± op/2 psi) 
for 4-specimen set (ref. 3) 

Before we determine any values, we can solve for t = Nt^ 



(17) 


Assume that en ~ 0.05 percent. Solving equation (17) gives t = 2.98x10^ 
minutes, and for N identical cycles each of = 3 minutes at temperature 
T, N = 990. 

The sensitivity of the results to variations in temperature and coating 
pulloff strength is illustrated in table II. 


Thermal transient effects . During the transient, the radial temperature 
gradient is strong and must be included. For constant heat flux the strain 
becomes (ref. 2) 
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o AaC|t 

e q = E = 2k ( 1 - v) 


( 18 ) 


which is effective primarily in the initial phases of the heating and cooling 
portions of the cycle. During the latter part of either cycle, as equili- 
brium is approached, the isothermal case would apply. (This would not be 
the case for an internally cooled specimen where a radial gradient would 
always be present.) Thus one can consider the heat flux to decrease. The 
allowable radial stress is 


rr 


E Aaqt 
Vr " 2k ( 1 - v)r 


(19) 


For q = 550 kW/m 2 , k = 0.069 cm-kW/m 2 -K, v = 0.25, E = 4.8xl0 4 MPa, 
and Ao = 6x10 _6 /K, the effective radial stress becomes 3.36 MPa and so 
the allowable pulloff stress, which we assume to be related to o rr , de- 
creases from 9 MPa (1300 psi) to 5.64 MPa (820 psi). This changes e to 
e + as follows. 



(20) 


The effects of the transient stresses in addition to the stresses created 
by creep, for variations in temperature and allowable stress, are illustrated 
in table III. 


Numerical Stress Model ing Codes 

In this section we present figures which define the thermomechanical 
response of a multilayered cylinder. For each numerical model and code we 
include a brief writeup relating to element type, mesh refinement, boundary 
conditions (mechanical /thermal ) , and material properties and discuss some of 
the results. 

ADINA (mechanical FE model) (ref. 34) . The element type employed was the 
small deformation, linear, thermoelastic material model which accounts for 
temperature-dependent data. The element is eight-node, two dimensional, 
isoparametric, and axisymmetric and is illustrated as figure 5(a). The axial 
z coordinate was chosen at 0.05 cm (0.020 in.) to maintain an acceptable 
aspect ratio (AR), thereby minimizing unrealistic distortion of the elements. 

The mesh was refined to determine the sensitivity of the model to varia- 
tions in node spacing. The coarse- and fine-grid models are illustrated in 
figure 1(a). The coarse mesh consisted of 18 elements, 4 in the ceramic, 1 
for the bondcoat, and 14 in the substrate. The refined mesh consisted of 66 
elements: 8 in the ceramic, 4 in the bondcoat, and 54 in the substrate for 

a total of 333 nodes and 666 degrees of freedom. 

The thermal boundary conditions discussed in the section HEAT TRANSFER 
ANALYSIS are enforced. The mechanical conditions are illustrated in figure 
5(a).*** In all cases the surfaces are traction free with o rr and 


***For the experimental case, there is always a temperature and a pressure 
distribution as well as a mechanical reaction to loadings and most likely 
flutter. While these loads are usually considered small, they will contri- 
bute to coating failures. 
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o rz = 0 at r = Rq and displacements at r = 0, u r = 0 because of 
symmetry and at least a node where u z = 0. 

Material properties and geometry: The temperature-dependent mechanical 

properties of the ceramic, bondcoat, and substrate used in this model are 
given in table I. It should be noted that the effective composition of the 
ceramic can be altered or graded, thereby adding other layers. Likewise can 
the bondcoat be altered and strain isolation provided. The substrate can 
also be graded. Thus multilayers can be used, but we broadly classify these 

as ceramic, bondcoat, and substrate. 

Some tentative results and conclusions: The calculated axial, radial, 

and circumferential stresses generated in a ceramic-coated cylinder are given 
in figure 6. These stresses result from the applied thermal and mechanical 
loadings described in the preceding sections. In this case, we have confined 
our attention to the material interfaces (figures 1 and 6), which represent 
maximum changes in stresses ('kinks' in the strain locus), and to the first 
15 seconds of the heating cycle to minimize the information presented. 

(1) Runs involving both coarse and fine meshes (18 and 66 elements) yield 
essentially equivalent results, but proper selection of node spacing is 
paramount. 

(2) Initial transient mechanical results, heating, and cooling are highly 
sensitive to the description of thermal profiles and mechanical boundary 
conditions. 

(3) Severe discontinuities in the ceramic, bond, and substrate properties 
yield adverse stress gradients; one needs to explore more compatible choices 
of material properties and geometries to define 'best' material -geometry 
pairing. It is hoped that materials science will generate requisite material 
properties combining good life (fatigue, fracture, etc.) (refs. 35-37). 

(4) Although ceramic properties have a profound effect on the thermal 
profile, it has little mechanical influence on stresses in the substrate. 

It appears that ceramic stresses (for times > 5 sec) are 'going along for 
the ride' generated by substrate expansion. This would be especially true 
for a thin ceramic layer. 

MARC (mechanical /thermal FE model) (ref. 31) . In this model (figs. 1(b) 
and 5(b)) we have extended the planar axisymmetric model to include an axial 
portion of the cylinder. We have considered two types of boundary condi- 
tions, fixed and free, realizing that the heated section of the experimental 
cylinder, modeled herein, satisfies neither; however, the results will be 
bounded by these two classic boundary conditions. 

For the fixed end conditions, no axial displacement at z = 0 or z = 2L 
(free in r), the model simulates the infinite cylinder. For the free end 
conditions, fixed at z = 0 (free in r) with no constraints on displacement 
at z = 2L, the model simulates a short cylinder of length 4L symmetric about 
the plane z = 0. Some of these results are presented as figure 7 for fixed 
end planes and figures 8(a) and (b) for free end planes for the thermal 
loadings discussed in the section NUMERICAL THERMAL MODELING CODES. 

Figure 7 shows that during the heating phase, the radial stress changes 
from tensile, initially, to compressive near equilibrium. The axial stresses 
are compressive and large. The circumferential stresses are at first com- 
pressive but quickly change to large tensile values. (Recall that this is 
an elastic analysis without failures.) During the cooling phase the stress 
magnitudes decrease, but not as dramatically as with heating. 

Figure 9 illustrates how the radial stresses change with heating and 
cooling at four locations near the bondcoat-ceramic interface. As can be 
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seen, the maximum tensile stress moves radially toward the surface and can 
be thought of as a thermal-driven stress wave. Failure of a ceramic under 
these conditions would be expected to move from the interface to slightly 
into the ceramic. For thicker ceramics the failure plane should move even 
further toward the surface. 

Figure 10 illustrates the 'bending stress' applied to the now barrel- 
shaped model. The inside elements (fig. 1(b)) represent the average of 
stresses over that element, so the actual values at the nodal points would 
be higher at z = 0 and lower at z = L. The dramatic shift in the stresses 
from tensile to compressive results from the short-cylinder assumption and 
the applied thermal load (see the section NUMERICAL THERMAL MODELING CODES). 
The peak stresses occur within the ceramic, rather than at the interface. 
Experimentally, most failures appear to occur near the bondcoat-ceramic in- 
terface. If the model were homogeneous, subsurface cracking within the 
ceramic could be the dominant mechanism (refs. 39 and 40). 

To determine some effects of finite element aspect ratio (AR), the ele- 
ment length was increased to unity. For the fixed constraint case the re- 
sults do not differ to any significant degree. 

FEATS (ref. 39) . For this code the geometric conditions are the same as 
in the thermal modeling section, (figs. 1(c) and 6(c)), and the materials 
used were both superalloy and the 304 stainless steel. The code was run 
with fixed and free boundary constraints where the aspect ratio was kept 

0. 5 < AR < 3, i.e., az = 0.0025 cm (0.001 in.). These constraints should 
bound the calculated stress for the experimental cylinder - if it remains 
elastic. In addition, the aspect ratio was made very large, 40 < AR < 2000, 

1. e., z = 2.54 cm (1 in.). This tends to artifically stiffen the substrate 
yet enable the ceramic and interfaces to respond to large temperature gradi- 
ents. It also artifically places the end constraint somewhere between fixed 
and free. Although the procedure is certainly questionable, the results 
correctly approach those theoretically predicted for homogeneous cylinders 
(ref. 40). 

The temperature profiles were generated by SINDA using the thermal load- 
ing discussed in the section NUMERICAL THERMAL MODELING CODES. The numerical 
stress values were generated by FEATS. As can be seen, the radial stresses 
are nearly the same for all three constraints, while the circumferen- 
tial stresses appear bounded by the free and fixed end constraint 

values. The axial stresses o zz follow these same trends. An exception 
occurs at the interface, where the fixed case stresses are large and com- 
pressive and the artificial case stresses are tensile. 

In general, there is a much larger variation in the stress profiles pre- 
dicted by ADINA, MARC, and FEATS than in the predicted thermal profiles, 
primarily because stresses are more sensitive to modeling methods than is 
the thermal analysis. 

SUMMARY 

Various levels of sophistication in both the modeling and the codes have 
been considered. The modeling ranged from one-dimensional constant proper- 
ties and first-order difference schemes with constant heat transfer coeffi- 
cients to temperature-dependent heat transfer and material properties used 
with an eight-node (serendipity) finite element method. Five numerical codes 
were involved in producing three sets of thermomechanical behaviors: ADINAT- 

ADINA, MARC, and SINDA-FEATS. The results of these numerical codes and a 
simplified analysis are as follows. 
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1. For the heating and cooling of a multilayered cylindrical heat sink 
specimen, it is necessary to include temperature-dependent fluid and material 
properties and coupling of the heat transfer from the hot gas stream to the 
specimen. 

2. Large thermal gradients, such as occur in multilayer-coated specimens, 
effect large thermoelastic strains, which for the elastic case give rise to 
high component stresses. 

3. During the heating or cooling transient, large thermal gradients exist 
in the skin of the coating and enhance the probability of failure through 
crack propagation. 

4. At steady operating temperatures, a simplified analysis indicates 
that inelastic behavior, creep, could play a significant role in specimen 
coating failure. It is suggested that the analysis be verified experimen- 
tally and through the use of numerical modeling of the multi layered-coated 
specimen. 

5. The temperature predicted using one-dimensional finite elements and 
finite differences has a variation of about 60 K. At this time, such varia- 
tions are thought to be the result of thermophysical properties and modeling 
procedures. It is also possible that there may be some variations due to 
solution techniques. 

6 . Coating failures usually occur at the ceramic-bondcoat interface. 

The calculated stresses at this interface, amplified by differences in 
thermal gradients, possess similar trends but the magnitudes differ. In 
general the elastic radial stresses are significantly less than the adhesive/ 
cohesive strength (at room temperature); however, consideration of inelastic 
behavior may enhance the calculated radial stresses. These problems are yet 
to be resolved. 

7. A major difficulty lies in characterizing problems which are neither 
plane strain nor plane stress. The desired result appears bounded between 
those for the fixed and free constrained disk-shaped surfaces of the model. 

An artificial constraint seems to provide results between these classic 
constraints but has not been explored in any detail. In general the 
application to the thermal problem is in violation of the field equations; 
nevertheless it is used. 

8 . The long- and short-cylinder behavior implies that the geometry 
either expands uniformly in the radial direction or it tends to 'barrel' at 
the plane of symmetry. Such behavior is characteristic of bending under 
thermal loads, but neither situation represents the finite cylinder. The 
next phase should look at the finite cylinder including axial and perhaps 
circumferential temperature gradients. 
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TABLE I. - REPRESENTATIVE THERMOPHYSICAL PROPERTIES 


Property 

Temper- 

ature, 

K 

Bondcoat 

(NiCrAlY) 

Ceramic coat 
(Zr02-Y 2 0 3 ) 

Substrate 

304- 

stainless 

Thermal 

283 

0.05769 

0.00721 

0.13848 

conduc- 

tivity. 

1352 

0.16589 

0.02164 

0.3346 

W/cm-K 

1770 

0.2082 

0.02728 

0.41129 

Density, 

kg/crrw 

All 

0.00698 

0.00567 

0.00782 

Specific 

283 

490.5 

490.5 

406.1 

heat. 

1367 

711.8 

670 

544.3 

J/kg-K 

1770 

794 

736.7 

595.6 

Modulus of 

283 

12.4xl0 4 

4.48xl0 4 

19.58xl0 4 

elasticity. 

811 

11.72 

3.21 

15.65 

MPa 

1144 

4.42 

2.41 

11.65 


1367 

a 0.0 

1.88 

8.55 

Thermal 

283 

12.6x10-6 

7.86xl0 -6 

17.1x10-6 

expansion. 

811 

13.86 

10.08 

18.59 

1/K 

1144 

16.2 

11.52 

19.53 

1367 

18 

12.24 

20.16 


( a Probably clxlO^) . 


17 










TABLE II. - ESTIMATE FOR NUMBER OF CYCLES TO FAILURE FOR VARIATIONS IN 
TEMPERATURE AND ALLOWABLE STRESS AT TWO DIFFERENT HEATING TIMES 


A1 lowable 

Number of 3- 

-minute cycles to failure 

Number of 57-minute cycles to failure 

stress 

Calculated at inter- 
face temperature- 

Experimental 
24-specimen 
data set 
(ref. 7) 

Calculated at inter- 
face temperature- 

Experimental 
4-specimen 
data set 
(ref. 2) 


1283 K 
(1850° F) 

1311 K 
(1900° F) 

1283 K 
(1850° F) 

1311 K 
(1900° F) 

°p ■ ”p - s o 

1.3xl0 4 

200 

810 

315 

5 

50 

o = o 
P P 

6xl0 4 

990 

1330 

1350 

22 

180 

a = o + S n 
P P o 

2xl0 5 

3340 

1870 

3585 

58 

250 
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TABLE III. - ESTIMATE FOR NUMBER OF CYCLES TO FAILURE FOR VARIATIONS IN 
TEMPERATURE, AND ALLOWABLE STRESS AT TWO DIFFERENT HEATING TIMES 
INCLUDING EFFECTS OF TRANSIENT HEATING 


A1 lowable 

Number of 3-minute cycles to failure 

Number of 57-minute cycles to failure 

stress 

Calculated at inter- 
face temperature- 

Experimental 
24-specimen 
data set 
(ref. 7) 

Calculated at inter- 
face temperature- 

Experimental 
4-specimen 
data set 
(ref. 2) 


1283 K 
(1850° F) 

1311 K 
(1900° F) 

1283 K 
(1850° F) 

1311 K 
(1900° F) 

°p - ”p - s 0 

920 

15 

810 

— 

— 

50 

0 = 0 
P P 

2.1xl0 4 

335 

1330 

265 

4 

180 

a =o + S n 

p p 0 

8.6xl0 4 

1375 

1870 

1230 

20 

250 
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Figure 1. - Finite element idealization of a multilayered cylindrical rod, thermal loading. 
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(a) Experimental (b) Model. 

(assumed). 

Figure 4. - Sketch of the assumed stress -time relations for a 
ceramic coated cylinder exposed to a 0. 3 mach burner flame. 
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C. SKETCH OF THE 105 ELEMENT MODEL USED IN FEATS, WITH THE 
FIXED, FREE AND ARTIFICIAL CONSTRAINTS 


Figure 5. - Finite element idealization of a multilayered cylindrical rod, mechanical loading. 
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Figure 6. - Stress distributions In a multilayered cylindrical rod subjected to the thermal and mechanical loadings of sections 3 and 4, ADINA. 
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Figure 9. - Radial stress distribution for elements adjacent to the 
plane of symmetry (inside elements) and those on the surface 
(outside elements). See Fig. 1 (b) for element numbers. 
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Figure 10. - Stress-time distributions for free, fixed, and artifi- 
cial end constraints. 
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